Long-range surface plasmon polariton excitation at the quantum level 
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We provide the quantum mechanical description of the excitation of long-range surface plasmon polaritons 
(LRSPPs) on thin metallic strips. The excitation process consists of an attenuated-reflection setup, where effi- 
cient photon-to-LRSPP wavepacket-transfer is shown to be achievable. For calculating the coupling, we derive 
the first quantization of LRSPPs in the polaritonic regime. We study quantum statistics during propagation and 
characterize the performance of photon-to-LRSPP quantum state transfer for single-photons, photon-number 
states and photonic coherent superposition states. 
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I. INTRODUCTION 

Plasmonics [ 1] is a rapidly growing area of research based 
at the nanoscale that is currently experiencing intensive stud- 
ies by researchers from many areas of the physical sci- 
ences yfl. Plasmonic -based nanophotonic devices using sur- 
face plasmon polaritons (SPPs) have recently started to attract 
much interest from the quantum optics community for their 
use in quantum information processing (QIP) 0, B la] • At 
present, it is essential that practical techniques are properly 
developed for efficiently generating and controlling plasmonic 
excitations at the quantum level. In order to do this, a rigor- 
ous quantum mechanical model must be included for describ- 
ing how photons and different forms of SPPs interact. With 
a clear description and theoretical understanding of these in- 
teractions, the rapid development of novel QIP applications, 
using nanostructured devices based on linear and nonlinear 
plasmonic effects will become possible. 

In this work, we adapt and extend techniques recently in- 
troduced by us [8] to provide the first quantum mechanical 
description of the coupling between single-photons and SPPs 
on thin metallic strips, also known as long-range SPPs [9] 
(LRSPPs). Here, an attenuated-reflection (ATR) setup is de- 
scribed, that has so far only been considered for classical LR- 
SPP generation [10]. In order to introduce the Hamiltonian 
for the interaction, we derive the first quantized description of 
the LRSPP fields in the polaritonic regime U. We find that 
high quantum efficiencies can in fact be reached for photon- 
to-LRSPP wavepacket-transfer upon appropriate modification 
of the ATR geometry. We comment on the extent to which 
the excited LRSPPs preserve quantum statistical features of 
the original photons as they propagate along realistic metal- 
lic strips. We then characterize the performance of photon- 
to-LRSPP quantum state transfer, focusing on an informative 
example of the transfer of coherent superposition states Hill . 
Recently, we have become aware of an experimental effort to 
transfer a similar type of nonclassical field into a LRSPP d. 

The benefits of exciting LRSPPs in this configuration 
compared to the previously studied standard single-interface 
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SPPs H [Ht] include a significant increase in the propagation 
length, together with the support for both transverse mag- 
netic (TM) and transverse electric (TE) polarization degrees 
of freedom, given the correct lateral width and thickness of 
the metallic strip Hill . Therefore they have the potential to 
open up a wider variety of QIP applications, where this type 
of additional flexibility is necessary. The work we present 
here provides a valuable description of the physics of photon- 
SPP coupling in multilayer geometries at the quantum level 
and the new methods we have developed specifically for this 
task should be well-suited to other complex SPP excitation 
scenarios, such as grating [15] and end-fire Hill techniques. 

The paper is structured as follows: In Section II, we provide 
the quantum mechanical description of LRSPPs. We also in- 
troduce the ATR setup used for the excitation of LRSPPs with 
single photons and the coupling Hamiltonian for the fields. In 
Section III we analyze this coupling for single modes of the 
system as well as for wavepackets involving single-photons 
and photon-number states. From this analysis we determine 
the transfer efficiencies of photons to LRSPPs over a range of 
input frequencies. In Section IV we then examine the extent 
to which quantum statistics of the injected photons are pre- 
served during transfer and propagation of the excited LRSPPs. 
In Section V we characterize the performance of photon-to- 
LRSPP quantum state transfer, providing an illuminating ex- 
ample of the transfer of coherent superposition states. Finally, 
Section VI summarizes our main results. 



II. EXCITATION SETUP 

LRSPPs are nonradiative electromagnetic excitations asso- 
ciated with electron charge density waves propagating along 
the interfaces of a dielectric-metal-dielectric configuration [9[. 
In Figs. Q] (a) and (b) we show the ATR setup and geometry 
utilized for single-photon excitation of LRSPPs. It consists 
of four layers: a metallic layer (with permittivity £3 = e m ), 
two dielectric layers (with €2 = £4 = 1), and a prism (ei). 
We consider the metal as silver in this work only to illus- 
trate our main results, with the theory developed supporting 
a more general setting. For LRSPP excitations, due to the col- 
lective nature of the electron charge density waves, a macro- 
scopic picture of the electromagnetic field produced is appro- 
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priate J^,[T3l. Upon quantization, LRSPPs are found to corre- 
spond to bosonic modes. A brief outline of this quantization 
is given in the Appendix. It is well-known that classically 
two types of TM surface modes are found for the geometry 
depicted in Fig. Q] (b): An antisymmetric mode, denoted by 
eigenfrequency lj + and a symmetric mode denoted by oj~. 
The quantized vector potential in the continuum limit for these 
modes propagating along an air-metal-air interface in the x di- 
rection, as shown in Fig. 1 (a), is found to be 

/>OC 

A$ PP (r,t) oc / do; ± (Ar ± (a; ± )W)- 1 / 2 x 
Jo 

[cp{r,uj)e^ luJ±t b(uj ± ) + h.c}. (1) 

Here A^ ± (w ± ) is a frequency dependent normalization and 
W is the beam-width lfl8ll . For both co + and uj~ the S(w)'s 
(W(u])'s) correspond to bosonic annihilation (creation) oper- 
ators which obey the quantum mechanical commutation re- 
lations [b(u>),b'(u>')] — S(u> — u)'). The modefunctions are 
given by 

± (r,w ± ) = e ik ' r [(k- {ik/u Q )i)e v ° z d{~z) 

+ (1 - v m /e m Vo)[Qs.+ (ik/is m )z)e~" mZ 
T(k- {ik/ Vm )i)e^ z ' d ^]d{z)d{d 1 - z) 
T(k + (ik/is )z)e~^ z - d ^tf(z - di)], (2) 

where the wavevector k = fcx, d(z) is the Heaviside step 
function and the decays into the metal and air are parameter- 
ized by v 2 n = k 2 - (a;*) 2 e m /c 2 and Uq = k 2 — (w ± ) 2 /c 2 
respectively. The dispersion relation between ui and k is 

e -v m d x _ ±(y m + e m i/ )/(i/ m - e m u ), (3) 

where d\ is the thickness of the metallic strip. The solutions 
of this equation, given by ui (k), correspond to two different 
types of coupled plasma excitations at the metal-dielectric in- 
terfaces 2/3 and 3/4 shown in Fig.[T](b), which oscillate syn- 
chronously out-of-phase (+) and in-phase (— ). The depen- 
dence of these solutions on the wavevector and slab thickness 
are shown in Fig. |2] (a), where silver has been chosen as an 
example having permittivity 11201 e m (w) = 1 — oj 2 /ui 2 + Se r , 




FIG. 1 : (Color online) Single-photon excitation of LRSPPs using 
attenuated-reflection, (a): A photon wavepacket is injected into the 
system at a specific angle 8, with a prism mediating an interaction 
between the photon and LRSPP modes. The minimum prism size is 
diffraction-limited, (b): The ATR excitation geometry considered, 
(c): Transfer process for the photon and LRSPP mode operators. 
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FIG. 2: (Color online) (a): Dispersion relations for the to excita- 
tions as a function of k and metal thickness di . The inset shows the 
curves cu + (k) (dashed) and u>~(k) (solid) corresponding to d\ = 20 
nm, and the same curves merged at d — 100 nm (dashed-dotted). 
(b): Coupling angle 9 as a function of frequency and metal thick- 
ness di. The inset shows the coupling angle 6 for uj + {k) (dashed) 
and u)~(k) (solid) at d± — 20 nm, and the same curves merged at 
d — 100 nm (dashed-dotted). 

with uj p = 1.402 x 10 16 rad/s and 5e r = 29w 2 /w 2 . For a 
fixed value of d\, they evolve above (+) and beneath (— ) the 
dispersion relation known for SPPs at a simple air-metal in- 
terface |1], and move closer to this curve as either k or d\ 
grows |9y, approaching the limiting value uj sp (where e m = 
— 1) as k — * oo for any d\. While only the oj + excitations have 
long-range propagation lengths, due to damping (addressed 
in detail in Sec. IV), here we consider both excitations as 
'long-range' in order to give a more complete description of 
the physical system. 

Following the diagram depicted in Fig.[T](a), let us consider 
an incoming photon propagating in the air with direction given 
by the unit vector k' = sin (9x + cos 6z and corresponding 
wavevector k' = fc'k' = k' x x + k' z z. The vector potential is 
given by Q 

A P (r,t)(x / dco(ojA)-^ 2 [e ik '^'- T h- iwt a(uj)+h.c]. 
Jo 

Here, A is the beam-width, with the a(w)'s and a'(u])'s sat- 
isfying bosonic commutation relations. The impossibility to 
fulfill the mode-matching conditions between the branches 
cj ± (/c) and the incoming photon beam, with dispersion rela- 
tion Lu(k') — ck' x / smQ (shaded region of inset in Fig.|2](a)), 
can be overcome in the ATR configuration by using a prism 
with dielectric constant ei > e 2 = 1, placed at a distance 
di over the surface of the metal. This modifies the dis- 
persion relation of the incoming photon beam to uj(k') = 
ck' x /(y/ei sin 0) (solid line of inset in Fig. |2]a). For 9 greater 
than the critical angle, an evanescent photon field is created 
below the prism surface due to total internal reflection lfl3ll . 
This provides a mechanism for achieving the coupling be- 
tween the photon and LRSPP, as the x component of the 
transmitted wavevector remains unchanged (see Fig.|2](b) for 
mode-matching 9 values). In order to model the ATR geome- 
try, we make use of the 4-layer (4L) configuration in Fig.Q](b). 
The modefunctions are given by 

*(r, uj) = rtp(r, uj)-d{-{z + d 2 )) + rip(r, u)d{z + da), (4) 

where r and r denote reflection and transmission coefficients, 
|r(w)| 2 + |r(u;)| 2 = l,Vw. These coefficients, together with 
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the modefunctions if>(r,oj) and tp(r,uj) are determined by 
solving Maxwell's equations for the incoming photon field, 
as in the classical coupled mode approach. In what follows, 
we will show how the r and r coefficients are combined with 
the quantum mechanical operators associated with the modes 
to derive the quantum coupling model. The modefunction 
■j/>(r, uj) possesses a real component of the wavevectorin the z 
direction and therefore cannot meet the mode-matching con- 
ditions required for coupling to LRSPPs. Only the modefunc- 
tion ip(r, u>) is involved in the coupling of photons to LRSPP's 
and we have 



7 TO (z-rfi) 



+(n 3 e 7mZ + K 4 e 
+« 5 e-^ (z - dl) tf(z-di)], 



(5) 
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where the Kj's are vector- valued functions related by bound- 
ary conditions at the interfaces, while 7^ = (k' x ) 2 — e m ui 2 /c 2 
and 7o 2 = (£4) 2 -„ 2 /c 2 . 

Within a linear response regime [8], the process of cou- 
pling between the photon field and the plasmon field can be 
described in the Heisenberg picture by a transformation ma- 
trix T(uj) as 



(6) 



with |a(w)| 2 + |/3(oj)| 2 = l,Vw. The transfer process is de- 
picted in Fig. Q] (c). The applicability of a linear approach 
is fully justified in this context, as we are interested in the 
description of the excitation of LRSPPs by a weak intensity 
photon field lfl9ll . The coefficients of T(uj) are determined 
through the overlap of system modefunctions, while the com- 
mutation relations of the quantum operators d(u>) and b(uj) 
properly define the structure of T(u>) as a valid unitary quan- 
tum transfer matrix Bill . The operators b in / out are associated 
with the in/out LRSPP modefunctions dr^ (r, uj) in Eq. @, i.e. 
bi n = b, whereas a in / out are associated with the in/out mod- 
efunctions \&(r, uj) of the 4L configuration. Here we assume 
that the photon field experiences negligible losses as it enters 
the prism and set di n = a. Thus, we have 



/3*H = -t(uj)5(uj - uj')S(k - k' x ) dz (M?(w)) 



X<^(z,u) • (jV 2 (uj')y 1/2 xP(z,uj') , 



-1/2 



(7) 



with ± (z, uj) and ip(z, uj') denoting the z~ dependent part of 
the functions in Eqs. (fJJ and (0, with normalization factors 
A/fk and N 2 respectively. Expression ^ can be obtained us- 
ing classical coupled mode theory. However, how the value 
of (3(uj) enters into the quantum coupling model of Eq. (O 
is determined by the commutation relations of the mode oper- 
ators a(uj) and b(uj) which define T(ui). A derivation of the 
coupling based solely on a classical electrodynamics approach 
is unsuitable when one considers properties that are explic- 
itly dependent on the quantum operators associated with the 
modes of the excitations. For instance, a description of the 
excitation of LRSPPs by an n-photon state of light would not 
be possible. This issue is discussed in more detail in Sect. IV. 



In describing the coupling process so far, we have used the 
Heisenberg picture. However, in many practical situations, it 
is more convenient to work in the Schrodinger picture, where 
the coupling is described by the following Hamiltonian 



H s = / duhwa^ (uj)a(uj) + / dujhujb 1 (u)b(uj) (8) 
Jo Jo 

/>oo 

+ih / du)[g(uj)a) (ui)b(ui) — g* (u>)P (u>)a(u>)], 
Jo 

with coupling coefficient JH g(uj) = e iar s/3(w) sin" 1 \/3(cu)\. 



III. PHOTON-LRSPP TRANSFER 

In this Section we investigate the optimization of the cou- 
pling coefficient g(uj) for a given range of parameters d\, d,2, 
and uj. For clarity we will handle the dependence on all three 
variables by first optimizing the coupling over d 2 for any pair 
of values of d\ and u>, and then optimizing over d\. Here, 
we use a prism with t\ = 1.51 and silver as an example, 
with the phenomenologically-derived i20ll dielectric function 
e m (u) = l-cj 2 /(cj(w + ir)) + <5e m ,wherer = 6.25xl0 13 
rad/s and Se m = 5e r m + iSe l m with 5e l m = 0.22. In deriv- 
ing the modefunctions ± (r, uj), we have assumed negligible 
damping in the metal for the surface plasmon field during the 
transfer process, allowing us to treat damping as the LRSPP 
propagates separately from the excitation process. 

For the optimization, some additional restrictions must be 
taken into account. First, under realistic conditions, the in- 
coming photon field does not constitute a genuine monochro- 
matic wave, but instead consists of a wavepacket with, for 
instance, a well-localized Gaussian distribution in frequency. 
This implies that both uj + and uj~ surface plasmons are sus- 
ceptible to being excited by the incoming wavepacket. For 
example, if the incidence angle 9 and the parameters, d\, d 2 , 
and uj, are set in order to achieve the excitation of one of 
the surface plasmons, either uo + or uu~ , it is possible that the 
other one might also be excited, due to the bandwidth of the 
wavepacket. To limit this effect we introduce a bandwidth pa- 
rameter, 
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Here, Aw is the bandwidth of an incoming Gaussian 
wavepacket, which we set as Auj = 3.02 x 10 13 rad/s, and 
is centered on uj — ui^ for £? ± . For a set d\ and uj ± , such 
that B ± > 1, the possibility to excite both surface plasmons 
simultaneously is very low. The values of B ± are shown in 
Figs. [3] (a) and (b) and indicate that only large values of d\ 
and low uj suffer from the possibility of simultaneous excita- 
tion. Note that the region corresponding to large frequencies 
and low d\ in Fig.[3](b) has been subtracted. This is because 
the dispersion for the symmetric excitation uj~(k) can never 
reach these frequencies for the range of d\ considered (see in- 
set of Fig. [2] (a)). This is applicable to all the plots for the 
symmetric excitation. 

A second restriction related to the optimization of the cou- 
pling parameter is the extent to which the LRSPP modefunc- 
tions ± (r,o;) penetrate into the prism. If the LRSPP field 
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FIG. 3: (Color online) Values of the parameters B , P , and C 
corresponding to the optimized coupling g ± (iv) over the separation 

d 2 . 



penetrates too much, the modefunctions should be modified 
to include the presence of the prism. In order to check the 
validity of using </> ± (r, to), we introduce a penetration factor, 
P ± = 2/v^d.2, which depends on the three parameters d%, d^, 
and lu. Here, P ± < 1 signifies that (/> ± (r,cj) at z = — d,2 is 
less than 2% its maximum value. The dependence of this fac- 
tor on d\ and to is depicted in Figs.[3](c) and (d). The values 
of P^ shown correspond to the optimized coupling coefficient 
g(u) (over the parameter c^) for any pair of values of d\ and 
uj. Since the value of g(uj) could deviate significantly from 
the true coupling in regions where P^ is large, due to the 
weak approximation of </) ± (r, u) to the true modefunctions, 
we must ensure g(uj) meets the condition P ± < 1. 

A final restriction concerns the coupled nature of the LR- 
SPP field. LRSPPs originate from the existence of coupled 
SPPs at both dielectric-metal interfaces of the metallic strip. 
For any finite value of the metal thickness, d%, it is always 
possible to find a solution of the dispersion relation Eq. (01. 
However, the interaction strength between both SPPs decays 
exponentially as d\ grows and the LRSPP evolves into two 
single SPPs. Due to a lack of symmetry during the single- 
photon excitation in the 4L configuration, the incoming field 
may be concentrated on the nearest metal surface to the prism, 
without reaching the other side. If the metal thickness is large 
enough, then it becomes impossible to excite any LRSPP. In 
order to account for this decoupling, we introduce a coupled- 
surfaces parameter, — 4/v^di. This factor quantifies the 
extent to which the field penetrates into the metal (see Eq. (ffj) 



to maintain a coupling between both SPPs, depending on the 
parameters d\ and u>. In Figs. [3] (e) and (f) the parameter C ± 
remains above 1 over the entire range considered, except for 
large values of d\ and low ui. Note that the restrictions and 
parameters introduced here also apply to the classical case, as 
they depend only the classical modefunction structure. 

While the penetration restriction was introduced previously 
in the context of single interface SPP excitation |8|], the band- 
width and coupling restrictions are new parameters emerging 
directly from this investigation of a multilayer configuration. 
With all three of these important restrictions properly iden- 
tified for the system, we are now in a position to correctly 
optimize the value of the coupling g(oj). In Figs. |4](a) and 
(b) we show the normalized coupling parameter, |g ± (tj)| = 
2\g ± (u)\/n, after being optimized over the variable c?2 and 
subject to the restrictions B ± > 1, P ± < 1, and C ± > 1. 
The regions where these restrictions could not be met have 
been subtracted. Here a value of \g \ = 1 (0) corresponds 
to the transfer of a photon field to a LRSPP with unit (zero) 
probability. The highlighted paths correspond to the maxi- 
mum coupling achievable. These paths are shown separately 
in Figs. [4] (c) and (d), plotted as a function of the frequency, 
together with the corresponding values of the parameters d\ 
and c?2- A maximum value of |<? + | = 0.9 is achieved for the 
uj + excitation, whereas the ui~ excitation shows a flatter be- 
havior, reaching \g~ \ = 0.8. It is interesting to note that as d\ 
increases, the optimized couplings in Figs. 0(a) and (b) move 
closer to those found for the photonic excitation of SPPs on 
a single interface 18], but with a remarkably lower efficiency. 
Under such conditions (c?i — > oo), it is only feasible to excite 
one of the SPPs and the excitation can no longer be consid- 
ered a single LRSPP (C ± — * 0). This transition from multi- 
layer coupled quantum excitation to single interface excitation 
should be an important factor to consider when designing op- 




FIG. 4: (Color online) Optimization of the coupling over the prism 
separation d 2 , subject to the constraints B > 1, P < 1, and 
C > 1. (a) u + excitation, (b) uj~ excitation. Panels (c) and 
(d) show the values of coupling parameter 1^*1 (dashed line), metal 
thickness d\ (solid), and prism separation d 2 (dot-dash), after numer- 
ical optimization: (c) lo + excitation, (d) u>~ excitation. For di, by 
increasing the resolution of the numerical calculation, one finds the 
surface plots in (a) and (b) become smoother and the points in (c) 
and (d) tend toward the best fit line. 
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timal quantum excitation methods in the context of multiple 
interfaces. 



IV. LRSPP PROPAGATION 

As excited LRSPPs propagate along the metal surface they 
experience loss due to finite conductivity of the metal and sur- 
face roughness. This results in heating and radiative losses 
respectively l23ll . For a reasonably smooth surface, thermal 
loss is the main source of damping |l|]. In order to include 
this loss mechanism for the LRSPP excitations, we follow a 
standard phenomenologic al ap proach, using a bath of quan- 
tized field modes 11221 124 12511 . The main advantages of this 
model are its simplicity and that it leads to the same physi- 
cal conclusions as a more rigorous derivation B26I1 . We con- 
sider an array of N = x/ Ax discrete, equally spaced beam- 
splitters, as depicted in Fig. (a). The upper ports repre- 
sent the spatial evolution of the operator for a propagating 
LRSPP, with input b out (oj) and output b^ ut (u>) after a dis- 
tance x. The lower ports consist of a bath of field exci- 
tations, di(uj), i — 1,...,N, which are independent and 
satisfy quantum mechanical bosonic commutation relations 
\c n {uj), & m {u)')\ — (>„ m (5(fjj — cj'). After applying successively 
the beamsplitter transformations, together with the continuum 
limit TV — > oo, Ax — > 0, such that c m (uj) — ^ \J Axc(uj, x') 
and 6 mn — > AxS(x — x 1 ), the operator of a damped LRSPP at 
point x can be written as ll22ll25ll 

b°Jd) = e iKx b 0Ut (u>) +iV2^ I* dx'e iK( - x - x ">c(u;,x'), 

Jo 

with [c(u), x), c>(ui', x')] = 5(u! — oj')5(x — x 1 ). Here we 
have introduced the complex wavenumber K = K(lo) = 
k(u))+in(Ld). This stems from solving the dispersion relations 
of the LRSPP excitations with the complex-valued dielectric 
function of the metal e m . Under the above conditions, the out- 
put LRSPP remains a bosonic excitation, with the appearance 
of the second term in b® ut (uj) from the bath of field excita- 
tions preserving this bosonic nature. We now define the time 
dependent creation and annihilation operators through the in- 
verse Fourier transform of the frequency dependent ones, for 
instance, b(t) = (27r) -1 / 2 / due~ iuJt b(uj). Although the lim- 
its of integration over the frequency are (—oo, oo), we are in- 
terested in the case of a narrow wavepacket centered on fre- 
quency uio with bandwidth Aui -C uJq, thus oj 6 (0,oo) can 
be taken. To calculate the mean flux of LRSPPs at a point x 
and time t after their excitation, we can write IT221 l25ll 

fouAt) = (b^ t (t)b° ut (t)) 

^ hj^j ^' (^ouMbout(Lo')) X 

e -(n(u)+K(u'))x e i[{k(u')-k{u))x-(u' -U))t] (o,) 

where we have used |H |H (&(uj,x)) = (c(u,x)) = 
((^(uj, x)c(lu', x')) = 0. Due to the small bandwidth of 
the wavepacket we are considering, the imaginary part of 
the LRSPP wavenumber remains essentially constant around 
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FIG. 5: (Color online) Phenomenological model to include the effect 
of damping for the LRSPP propagation 1231 . (a): Array of beam- 
splitters and bath of field excitations, (b) and (c): Normalized mean 
excitation number (m) = (m)/n as a function of the frequency and 
the distance traveled from the injection point. The optimal profiles 
obtained in Fig.|4]are used for the lo + excitation (b) and lu~ excita- 
tion (c). A detector efficiency of [i — 0.65 is used [8]. 

the central frequency, i.e. k(ui) w k(ujo) = kq, whereas 
the real part can be approximated by truncating its series up 
to first order, k(uj) = fc(oj ) + (uj — uj )v g ~ 1 (ujq) . Here 
Vq (wo) = dk(u>) I 'dui\u is the inverse of the group veloc- 
ity of the LRSPPs at the central frequency ui (for either u> + 
orw"). We then have 

(SLHWa/)) (10) 

= ^ 2KaX (bl ut (t R )b out (t R )) = e~ 2K °*f out (t R ). 

The mean flux of LRSPPs at point x and time t therefore 
equals that at x = and time t R = t — x/vg, but damped by 
a factor e~ 2K ° x due to losses incurred during the propagation. 
We now consider an incoming n-photon wavepacket state I22I1 
entering the prism with frequency profile £(w), given by 
\n 6 ) = (n!)- 1 / 2 ^)" |0). Here, a\ = fdcj£(u)tf(u) = 
J dt^(t)a^ (t) and for simplicity we use the Gaussian profile 
1(uj) = (27TCT 2 )" 1 / 4 exp[-i(w - uj )t Q - (u - uj ) 2 /Aa 2 }, 
where a — Aw/(2V2 In 2) and to is the time of injection. 
The mean number of LRSPPs that can be detected at point 
x is obtained by integrating f® ut (t) in Eq. ( TTOt over the time 
interval r = [to + x/vq — 3er, to + x/vq + 3<r], leading to 

(m) =JfZt(t)dt = ne~ 2 ^\PM\ 2 n, (11) 

where /i parameterizes the detector efficiency i27ll . and we 
have used the fact the LRSPP coupling does not vary appre- 
ciably over the bandwidth considered QSQ, i.e. (3{lo) ~ f3{wo). 
Figs.[5](b) and (c) show (m) = (m) /n, for oj + and uj~ respec- 
tively, along x for the optimized (3(lo) (and therefore g(co)) 
obtained from Figs. H](c) an d (A). It is interesting to see the 
positive effect that the structure of the lo + excitations has on 
damping reduction compared to cj~. 
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While the observables (m^) match well the behavior of 
their classical counterparts [1], the field intensity /, we must 
emphasize that the formalism presented here provides a more 
complete description of the LRSPPs; we are now able to in- 
vestigate the behavior of quantum statistics. In particular, in 
order to show that the LRSPP field has quantum character- 
istics, we need to consider the zero time-delay second-order 
quantum coherence function g( 2 )(0) at a fixed position l22tl . 
This observable is defined as (0) = ( : I 2 (t) :)/(: I(t) : ) 2 , 
where I is the intensity of the quantized field operator, : : 
denotes normal-ordering of the quantum operators and the 
expectation value is taken over the initial state of the field. 
It has been noted recently [8], that as the photon-to-surface 
plasmon transfer and propagation stages constitute an ar- 
ray of lossy beamsplitters 12711 . g^(0), which is equal to 
(m(m — 1))/(to) 2 for an incident n photon wavepacket, re- 
mains unaffected by the entire conversion process. This is 
to be expected [22], because at a beamsplitter with loss co- 
efficient 77 1 / 2 , the quantum observables (m) — > ry(m) and 
(m(m — 1)) — > r} 2 (m{m — 1)). Thus, the individual losses 
accumulated from the transfer and damping processes can- 
cel due to the form of g^(0). For a classical field 1 < 
5 (2) (0) < 00. Thus .g (2) (0) for a propagated LRSPP will al- 
ways lie in the classically forbidden region g^> (0) < 1. A 
Hanbury-Brown Twiss type experiment [28] could be used to 
measure g^ 2 \Q). Here, single-photon detection-based tech- 
niques could be employed. One might use an additional prism 
to convert the LRSPP excitation back into a photon and in- 
directly measure the signal with avalanche photodiode detec- 
tors. A more direct approach would be to use avalanche-type 
plasmonic detectors Il29ll embedded within the metal surface 
to directly probe the surface plasmon's excitation signal. In 
either approaches, the detection data from many identical ex- 
citation processes would be required in order to determine the 
overall quantum expectation value of g^ 2 \0). This repeti- 
tion technique is used frequently in quantum photonic experi- 
ments [ 30] and could be achieved easily by using a steady rate 
of single photons injected into the prism at set time intervals. 

V. QUANTUM STATE TRANSFER 

The ability to interconvert a quantum state between two 
kinds of physical system is an important requirement for QIP 
and is one of DiVincenzo's criteria [33] for quantum comput- 
ing and communication. As an example of efficient quantum 
state transfer from photons to LRSPPs using the theory devel- 
oped in the previous sections, we consider a superposition of 
coherent states [11]. These states reside in a single field mode 
consisting of superpositions of two coherent states of equal 
amplitude separated in phase by 180° and written as 

W)=N{\a)+e i «\-a)). (12) 

Here \±a) = exp(-|a| 2 /2) £^ =1 [(±a)"/v^!] \n), with 
the normalization TV = [2 + 2exp(— 2|a| 2 ) cosv?] -1 / 2 . When 
a is real, \ip) for tp = is orthogonal to the state with tp = tt, 
regardless of the size of a. Therefore one can use \tp) and 
\tjj) OT as a logical basis for QIP 03411 . For the remainder of 



the paper, we will focus on \ip) as our example basis state. 

Photonic superpositions of coherent states injected into the 
prism can be transferred to a quantum state of LRSPPs under 
the transformation matrix T(uS) given in Eq. ©. The unitary 
transformation defined by T(lo) acting on the input product 
state \^f) in = \ip) a . |0)b ire produces the state 

\®)out = -^(l" cos 3( w ))a ouf |-" sin 5(^)> bouf 

+ \-a cos g{u)) aout \asmg{u))) bout ), (13) 

where the coupling coefficient g(u) = sin |/3(a;)| is used, 
with phase $ absorbed into the definition of the incoming- 
outgoing fields 112 ill . Tracing out the unobserved photon mode 
a out of the photon-LRSPP system, the final state of the LR- 
SPP mode b out for a specific frequency to is given by 

Pb out — N' 2 (\a smg){a sin g\ + cq \a sin g) (— a sin g\ 

+cq |— a sing) (a sin g\ + \— a sin g) (— a sing|), 

(14) 

where cq = exp[— 2a 2 cos 2 g}. We now consider the damped 
propagation of these excited LRSPP coherent superposition 
states as they travel along the surface of the metal. This anal- 
ysis complements well the study performed in the previous 
section on amplitude damping of the LRSPP field. In contrast 
to before however, we now have a coherent superposition of 
an LRSPP excitation and we seek to characterize the effects of 
loss of coherence in this state, more commonly referred to as 
decoherence. Using the identity for the damped LRSPP field 
operator b„ ut in Section IV, and applying it to an excited co- 
herent superposition state wavepacket of central frequency ujo, 
as described by Eq. ( TBI , a straightforward calculation leads 
to 

Pb(x) = Af' 2 (\asmge- KoX )(asmge- K ° x \ 

+ cqc(x) \asing e~ K ° x )(-asmg e~ K ° x \ 
+ cqc(x) J —a sin g e~ K ° x ) (a sin g e ~ K ° x | 
+ |-asin 5 e~ Kf,x )(-asin ff e~ K ° x \) : (15) 

where c(x) — exp[—2a 2 (sin 2 g){l — e~ 2K ° x )]. It is eas- 
ily seen from Eq. (TT3T > that as LRSPP coherent superposition 
states propagate along the metal surface, the initial superposi- 
tion evolves into a statistical mixture of coherent states due to 
the factors cqc(x). As the coefficients of the off-diagonal el- 
ements of the density operator expressed in the coherent state 
basis vanish fastest, the initial mixture of the coherent state 
superposition tends toward a classical mixture (dephasing) at 
early times, eventually moving toward the vacuum state (am- 
plitude damping) at long times. The coefficients of the off- 
diagonal elements also indicate that the greater the value of 
a, the more quickly quantum coherences will decay through a 
dephasing type process of the LRSPP state 1 3711 . 

A characterization of the loss of coherence in a quan- 
tum state can be investigated using the von Neumann en- 
tropy [35] defined as Sy = — Tr(p In/3). The von Neu- 
mann entropy is a monotonic function of the linear entropy 
for a two-level quantum state of a single mode [36]. The 
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FIG. 6: (Color online) The entropy Sv as a function of frequency 
and distance traveled for fixed values of cat-state amplitude a. (a) 
and (c) correspond to the lu + LRSPP excitations with a — 2 and 
5 respectively, (b) and (d) correspond to the u~ excitations with 
a — 2 and 5 respectively. 



evaluation of the von Neumann entropy requires, in gen- 
eral, the diagonalization of the density operator p. Fortu- 
nately, the density operator of coherent superposition states 
with a E M can be decomposed into the orthonormal ba- 
sis |±) = A/±(|asin<7 e~ K ° x ) ± |— a sing e~ K " x )), with 
A/"± = (2±2e- 2Q2 ( sin2 ^ e_2 ' t0 ")- 1 / 2 . In this basis, Eq. (B) 
becomes 



VI. CONCLUDING REMARKS 

We have provided a fully quantum-mechanical description 
of the photonic excitation of LRSPPs using a versatile ATR 
geometry. In order to do this it was necessary to quantize the 
LRSPP field, which we included as an appendix. With this, we 
described the photon-LRSPP coupling mechanism by means 
of a linear Hamiltonian and optimized the coupling efficiency 
over a wide-range of parameters accessible to the setup. We 
found remarkably good transfer efficiencies. A phenomeno- 
logical model was then used to account for damping as the 
LRSPPs propagate, where the long-range behavior of the exci- 
tations manifested itself through quantum interactions with an 
environment. The effect of finite bandwidth for the incoming 
photon field on the coupling optimization and on the propaga- 
tion of a LRSPP wavepacket along the metal surface was also 
discussed. We studied the quantum statistics of the excited 
LRSPP fields and provided an outline of how one might exper- 
imentally investigate them. Finally, we characterized the per- 
formance of photon-to-LRSPP quantum state transfer, provid- 
ing an informative example of the transfer of coherent super- 
position states 111 ill . We found efficient transfer and analyzed 
the loss of decoherence in the states. The work presented here 
should be a useful starting point for future research into the 
practical design of novel long-range and multilayer plasmonic 
quantum-controlled devices based at the nanoscale. Applica- 
tions in this context include SPP-enhanced nonlinear photon 
interactions and SPP-assisted photonic quantum networking 
and processing. 



p b (x) = X + (x)p + (x) + X-(x)p-(x), 



(16) 



where X±(x) = 2 ^C a (1 ± cqc(x)) are eigenvalues cor- 
responding to eigenstates p±(x) = |±)(±|, together with 
A+(x) + X-(x) = 1. The von Neumann entropy is then sim- 
ply Sy(x) = — A+ln(A+) — A_ln(A_). In Fig. |6]we show 
the dependence of Sy on the frequency and the distance trav- 
eled for specific values of a. From Fig.|6]one can see that the 
greater the value of a, the more quickly the entropy increases 
toward unity, and thus the mixedness of the state increases, in- 
dicating a greater loss of coherence due to dephasing, an effect 
noted earlier ll37ll . Moreover, one can see from the left hand 
{u! + ) and right hand (cu~) columns of Fig. |6] that, for a given 
value of a, the entropy slowly increases for the u> + LRSPPs 
compared to the lu~ excitations as they propagate. This effect 
is related to the amplitude damping process observed in Sec. 
IV for the propagation of LRSPP-number states and is due to 
the smaller value of ko for a given mode frequency for the ui + 
excitations, a result of the structure of the modefunctions and 
their corresponding dispersion relation. According to Jeong et 
al. B38I1 mixed macroscopic superposition states can, in some 
cases, be more robust with respect to decoherence than their 
pure state counterparts. We expect this study to be useful in 
future work on the optimization of quantum state transfer of 
photons to LRSPPs and the consideration of different forms 
of SPPs on multiple interfaces. 
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APPENDIX: LRSPP QUANTIZATION 

Here we develop the quantization procedure of Elson and 
Ritchie [17] for the case of a thin metallic strip. An alternative 
approach based on the point-ion model for a dielectric slab has 
been used in Ref. 13911 to quantize the long-range plasmonic 
fields in the polaritonic regime. For quantization in the limit 
of short wavelengths, i.e. the non-retarded regime, the reader 
is kindly referred to Refs. JKH- 

Classical mode structure - We start with the geometry de- 
picted in Fig. 1 (b), where Maxwell's equations in terms of 
the vector potential A(r, i) lead to 



V 2 



c 2 dt 2 
e 



A = V - (V- A), 



V • A = -(n(r,t)-n (r)). 



(17) 
(18) 



Here e = e(r, oS) is a position and frequency-dependent di- 
electric function, n(r, t) is the electronic number density of 
the electron gas, no(r) is the static density in the undisturbed 
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electron gas and e is the electronic charge. We use the gauge 
4> = 0, where the electric field E = — A and magnetic field 
B = V x A. The classical energy residing in both the fields 
and electron gas is given by llT7ll 



n 



£0 
2 



1 + ?9 



eg 



2 (V x A) 2 



Here, § eg = §(z)§(d - z) 



(19) 

is a step function for the elec- 
tron gas located in the region < z < d, with t? de- 
noting the Heaviside function. In addition, u> p (r,t) — 
[n(r, t)e 2 /(eom*)] 1 / 2 is a position and time dependent 
plasma frequency, with m* the effective electron mass. Fol- 
lowing a linearized hydrodynamic approach lfl7ll32ll and tak- 
ing into account the location of the electron gas, the approx- 
imation n(r, t) « ^o( r ) = 0(z)i?(<2 — z)no is used. Corre- 
spondingly, we have e(r, 10) — $(—z) + $(z){}(d — z)e(u>) + 
d(z)'d(z — d), where e(u>) is a real- valued dielectric function 
for the metal of thickness d sandwiched by two layers of air 
with e = 1. Note that here we are considering an ideal case 
with no damping effects in the metal. This simplifies the quan- 
tization procedure. However, damping can be introduced at a 
later stage [8], as described in the main text. From Eqs. (fTTT i 
and (fT8l and the above considerations, we now have the clas- 
sical field equation 



dt 2 



0, 



(20) 



and V • A = in the region z 7^ {0, d}. Here, the usual 
conditions of continuity of the tangential components of the 
fields across the planes at z — and d respectively must be 
satisfied. To find the normal mode solutions for the system, 
we make the standard Ansatz 



A(r,i) =^A k (z)7V k (t) e l 



(21) 



where r = ix + yy is a vector parallel to the x — y plane, k = 
k x x + k y y and the associated eigenfrequency, denoted by w^, 
depends on k. The temporal amplitude iVk(i) is assumed to 
satisfy the oscillator equation, i.e. (d 2 /di 2 + ufyN^t) = 0, 
thus upon inserting Eq. (fJTJ into Eq. (|20l i, one obtains 



A±(z) = 0, 



dz 2 



i/ m Aj?(z) = 0, 



where — k 2 — w^e(wk)/c 2 and v 2 = k 2 - lo^/c 2 . 
Here the spatial amplitudes A k (z), AZ{z) and A™(z) cor- 
respond to fields in the z > d, z < and < z < 
d regions respectively. Consider the following solutions: 

! and AfO) = 



A+(z) = A+e-Mz-d), A-(z) = A, 



,u m {z-d) wheJ . e A + = a + k 



a k k + f3 k z and A™* = ag^k 



(3^ z. Here we 
focus on TM modes of the system due to boundary condi- 
tions ||9j,ll4j]. By requiring that the tangential components of 
E and B derived from A across the planes z — and d are 
continuous and V • A = elsewhere, one finds the oik's and 
/?k's are related to one another, with solutions existing only if 
e -vmd _ ±{y m + e(u!k)uo) / (v m — e(wk)fo) is satisfied. At a 



set thickness d there are two possible eigenfrequencies of this 
equation for a given k, which we denote Thus, there are 
two sets of coefficients for the a^'s and (3^s for a given k. 
The most general form of A(r, t) is then 



A± (r, t) = 4>t CO^k (ty k r + c.c, 



(22) 



where the eigenmodes are given by 

+ (1 - v m /e(u-k)vo)[(k + (ik/v m )i)e~ VmZ 
T(k - (ik/is m )iy^ z - d ^(z)d(d - z) 
T(k + {ik/u )2.)e- V0 ^-^d{z - d)], (23) 
and the time-dependent amplitudes are given by N^(t) = 



In the above, we have used ov as the free coef- 



ficient in the coupled boundary equations and absorbed it into 
the definition of 2Vjf . Due to the symmetry in the phases of 
the amplitudes in the 0jr (z)'s with respect to the center of the 
metal, the associated field modes are commonly referred [ 1, 9] 
to as antisymmetric (w k ) an d symmetric modes (w k ). 

Discretization and quantization - We now discretize the 
classical system and quantize it. The components of k are 
taken to be k x = 2ir£ x /L and k y = 2tt£ v /L, where i x and £ y 
are integers, with e satisfying boundary conditions at the 
planes x = ±L/2 and y = ±L/2. Substituting Eq. d22l) into 
Eq. ( fT9l one finds the total energy of the discretized classi- 
cal modes given by W± = £ k £oL 2 (uj± ) 2 Af± {N± N± * + 
N^*N^), where A/" k is a coefficient with dimensions of 
length. Using the correspondence with a quantized harmonic 

^) 1/2 k.± and 



oscillator MM, i-e. N, 



T± 



(H/2e L 2 uj^jv k ) -u k . 



±* 



(ft/2e £ 2 o; k t A/ k t ) 1 / 2 Sj c we have the Hamiltonian 



= -hu,±(b kt± bl ± +bl ± b K± ), 



(24) 



along with the vector potential converted to the operator 

A±(r,*) = 5](^/2e o i 2 a;±AA k ± ) 1 / 2 [0±(r)e-^ ± *6k, ± +^.c]. 



k 



where the creation 



Here we have </> k (r) 
and annihilation operators for the quantum excitations satisfy 

[&k,±,& k ',±] = ^k,k'- 

Continuum limit and beam-width - We now take the con- 
tinuum limit using the transformations ^ k — > (L/2iv) 2 f dk 
and 6k, ± — > (27r/L)6±(k), leading to 

27 



A ± (r,t) 



Ak{h/2e Q uj ± {k)N ± {k)) 1/2 x 



[^(r.kje" 



-iw ± (fc)t 



b±(k) + h.c}. (25) 



Next, the excitations propagating in the x direction have a 
beam- width W imposed in the y direction lfl8ll using J dk — > 
(2tt/W) J2k I dk x and setting k y = 0, with (5 2 (k - k') -> 

(W/2-rr)8{k - k') and S±(k) -> (W 1 / 2 /27r)6±(fc), giving 
A ± (r,i) = i- j dk(h/2e Wuj ± (k)M ± (k)) 1 / 2 x 

[^(r, k)e- lbJ±{k)t b±{k) + h.c}. (26) 
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Finally, we convert to the frequency domain using dfc — + 
[« G (w ± )]- 1 dw ± and b±(k) -> [u G (w ± )] 1 / 2 S ± (w ± ), where 
«g( w± ) = duj^/dk is the group velocity. This gives the 
quantized vector potential for the uj^ field as 

A ± (r,t) = i /da; ± (fi/2eoWa; ± u G (a; ± ) J V" ± (a; ± )) 1 / 2 x 



2tt . 



[0 ± (r,w ± )e-^ *6±(a; ± )+/i.c]. (27) 
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